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OO ■ We propose a scheme for determining a generalised scaling function, namely the 

Sudakov factor in a peculiar double scaling limit for high spin and large twist 

f— >) ■ operators belonging to the sl{2) sector of planar AA = 4 SYM. In particular, we 

OO , perform explicitly the all-order computation at strong 't Hooft coupling regarding 

the first (contribution to the) generalised scaling function. Moreover, we compare 
our asymptotic results with the numerical solutions finding a very good agree- 

kS ' ment and evaluate numerically the non-asymptotic contributions. Eventually, we 

C^ ' illustrate the agreement and prediction on the string side. 
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1 On the set-up 

The s/(2) sector of A/" = 4 SYM contains local composite operators of the form 

Tr{V'Z^) + ...., (1.1) 

where V is the (symmetrised, traceless) covariant derivative acting in all possible ways 
on the L bosonic fields Z. The spin of these operators is s and L is the so-called 'twist'. 
Moreover, this sector would be described - via the AdS/CFT correspondence [1] - by 
string states on the AdSs xS^ spacetime with AdSs and S^ charges s and L, respectively. 
Moreover, as far as the one loop is concerned, the Bethe Ansatz problem is equivalent 
to that of twist operators in QCD [21 [3]. u 

Proper superpositions of operators (11.11) have definite anomalous dimension A de- 
pending on L, s and the 't Hooft coupling A = Sir'^g'^: 

A = L + s + 7((7,s,L), (1.2) 

where '-/{g, s, L) is the anomalous part. A great boost in the evaluation of '-/{g, s, L) in 
another sector has come from the discovery of integrability for the the purely bosonic 
so(6) operators at one loop [1]. Later on, this fact has been extended to the whole 
theory and at all loops in the sense that, for instance, any operator of the form (II. ip is 
associated to one solution of some Bethe Ansatz-like equations and then any anomalous 
dimension is expressed in terms of one solution [5] . Nevertheless, along with this host of 
new results an important limitation emerged as a by-product of the on-shell (S-matrix) 
Bethe Ansatz: when the interaction range had been greater than the chain length, 
then unpredicted wrapping effects would have modified the same. More precisely, the 
anomalous dimension is in general correct up to the L — 1 loop in the convergent, 
perturbative expansion, i.e. up to the order g"^^'"^. Which in particular implies, - 
fortunately for us, - that the asymptotic Bethe Ansatz gives the right result whenever 
the L —>■ oo limit is applied. 

In this limit an important double scaling may be considered on both sides of the 
correspondence □: 

s-^oo, L^oo, j = - — = fixed . (1.3) 

ms 

In fact, the relevance of this logarithmic scaling for the anomalous dimension has been 

pointed out in [7j within the one-loop SYM theory and then in [6] and [8] within 

the string theory (strong coupling g ^ 1). Moreover, by describing the anomalous 



^A deep reason for that may be that one loop QCD still shows up conformal invariance, albeit 
integrability does not seem to hold in full generality (for instance, it apparently imposes aligned 
partonic helicities). 

^Actually, in string theory (semi-classical calculations) the A -^ +oo limit needs consideration 
before any other (cf. for instance [B] and references therein), thus implying a different limit order with 
respect to our gauge theory approach (cf. below for more details). 



dimension through a non-hnear integral equation [9] (hke in other integrable theories 
[To]), it has been recently conjectured the Sudakov scaling [8] 



^ig,s,L) = f{g,j)\ns + Oii\nsm, (1.4) 

thus indeed generalising to all loops the analogous result of [7]. |j Actually, this state- 
ment was argued by computing iteratively the solution of some integral equations and 
then, thereof, the generalised scaling function, f{g,j) at the first orders in j and g^: 
more precisely the first orders in g"^ have been computed for the first generalised scaling 
functions fn{g), forming 



oo 



f{9,j) = }^fn{9)r. (1.5) 

n=0 

As a by-product, the reasonable conjecture has been put forward that the two-variable 
function f{g,j) should be analytic (in g for fixed j and in j for fixed g). In [H] similar 
results have been derived for what concerns the contribution beyond the leading scaling 
function f{g) = fo{g), but with a modification (see also [12] for the general idea) which 
has allowed us to neglect the non-linearity for finite L and to end-up with one linear 
integral equation. The latter does not differ from the BES one (which cover the case 
j = 0, cf. the last of [5]), but for the inhomogeneous term, which is an integral on the 
one loop root density. In this respect, we have thought interesting the analysis of the 
next-to-leading-order (nlo) term - still coming, for finite L, from an asymptotic Bethe 
Ansatz -, as the leading order f{g) has been conjectured to be independent of L or 
universal (penultimate reference of [5], after the one loop proof by [7]). 

For instance, the linearity in L of this nlo term and other features have furnished us 
the inspiration for the present calculations. In this respect, we deem useful to briefly 
introduce in the next Section a suitable modification of the aforementioned method 
such that it applies to the regime (II. 3p (still for any g and j). In principle, we might 
in this way determine the generalised scaling function f{g,j), and also its constituents 
fn{g)- As a consequence of the full range of g we can work out not only the weak 
coupling expansion, but also the exact strong coupling evaluations (here only for the 
simplest case concerning fi{g)), and to perform an all g numerical work. Eventually, 
we shall highlight the concise style of the present publication, in that a next one with 
much more details shall be coming out rather soon. 

2 Computing the generalised scaling function 

Because of integr ability in A/" = 4 SYM, we have been finding useful [13] to rewrite 
the Bethe equations as non-linear integral equations, whose integration range is, even 
for massive excited states (holes and complex roots), (— oo, -|-oo) [TD]. More recently 



■^0((ln s) °°) means a remainder which goes faster that any power of hi s: Um (In s)'°0((ln s) 

s — >oo 

0,Vfc>0. 



[TT1[T2], we have developed a new technique in order to cope with an integration interval 
that may be kept finite, having in mind the peculiar case of the s/(2) sector with two 
holes external to all the roots [3 El [9]. With reference in the entire Section to [TT] . 
the non-linear integral equation for the s/(2) sector involves two functions F[u) and 
G{u,v) satisfying linear equations. Splitting F{u) into its one-loop and higher loop 
contributions, Fq[u) and F^{u) respectively, the latter has been shown to satisfy (cf. 
(4.10) and (4.11) of ^): 
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which holds when s ^ oo once introduced in equation (4.10) of [H] the Bethe roots 
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Here the hole contributions given by Sh{v) = {L — 2)6{v) + 0(L^/lns) and rih{v) = 
0{L^/ In s) depend on the logarithmic scaling to zero of the internal holes. Moreover, the 
solution of (12. ip yields the anomalous dimension in the large spin limit while neglecting 
the non-linear convolution terniQ 
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(2.3) 



where the last equality is a generalisation of the analogue in |14j . 

As we are just neglecting terms smaller than any inverse logarithm, 0((lns)~°°), 
we can efficiently consider the double scaling limit (II. 3p . Looking at (12. ip . the r.h.s. 
is made up of seven terms, among which the first is explicit and the last uninfiuential. 
On the contrary, the third and the fifth integral refiect two different kinds of important 
hole contributions. Nevertheless, thanks to the quoted scaling 6h{v) = {L — 2)6{v) + 
0(L^/lns) and T]h{v) = 0(L^/lns), these integrals enjoy an expansion about j = 



'''The actual positive densities would also have a factor —1/{2tt) in front [TT]. 

^The dependence bo{s) ~ s/2 + 0{s'^) comes out from the normalisation of the density <Jo{u) [TT]. 



^We also have some numerical evidence supporting this fact. 



(similar to the following one) contributing only with the Dirac delta function to fi{g), 
and the remainder (in particular the all f]h{v)) to the higher /„((?), n = 2, 3, ... . And 
last but not least, the second term determines the inhomogeneous (or forcing) term as 
a series about j = via the density cro(^); by virtue of the logarithmic expansion [3, |9] 
of the one- loop theory [11] : 
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In fact, this series shall follow from the solution (in the scaling (II. 3p ) of the linear integral 
equation for the one loop density o"o(u) (equation (4.9) of [H]), which determines all 
the functions 0n(w). As a consequence, the solution of ( 12. ip inherits the same form of 
the forcing term 
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Eventually, all the generalised scaling functions f{g,j) and fnig) may be computed, in 
principle, via the equality (12.31) . From now on, we will restrict ourselves to the first 
order beyond the j = theory [5], i.e. (7^ (u) and fi{g). 

3 The first generalised scaling function 

In [11] we have found a way to compute the Stieltjes integrals with measure induced 
by the one loop-density o"o(f) on the support {—bo, bo) via an effective density ctq (f) 
on support (— oo, +oo) 
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where the Fourier transform of the function ctq ° (w) reads explicitly 



(bo), 



a, 



(bo) 



(k) 



-Att 



e 2 cos -^ 



2sinh 



|fc| 



27r(L-2) 



_ifcl 
e 2 



2sinh 



|fc| 



(3.1) 



4:7iS{k)\n2, (3.2) 



and q{v) is a suitable test-function. Therefore, we can compute exactly the functions 
(j)o{u) and 0i(m) and thus write down the linear integral equations for (J^ (u) and cr^j (u) 
respectively. Of course, the equation for o"]^ (m) is the BES equation [5], whilst that 
about ajj (u) is a novelty. We may find still convenient to manipulate the Fourier 
transforms defining the even function 



^w^^^S'w 
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which has the important property (induced by (12. 3p ) 
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If we introduce for convenience the functions 

«r(^) = / —Jr{V2gh) ^, ar{g)= —Jr{V2gh) ^ 

and re-cast the integration inside the domain k>Q, the integral equation for s{k) takes 
the form 
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As anticipated, this equation shares the same kernel with BES equation [5] , but different 
forcing term. Its solution may be expanded in a series involving Bessel functions: 
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Now we want to derive linear equations for the coefficients s„. Upon expressing the 
scattering factor coefficients via 

dh Jn{\/2gh)J„,{\/2gh) 
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equation (13. 6p decomposes into the infinite dimensional linear system 

(OO CO 

-a2p{g) - 2 ^ Z2p,2m{g)S2m{g) + 2 ^ Z2p,2m-l{g)s2m-l{g) 
m=l m=l 

(3.10) 

/ s OO OO 

2^"^-^ = -a2p-l{g) - 2 ^ Z2p-l,2m{g)S2m{g) - 2 5]] Z2p-l,2m-l{g)s2m-l{g) , 

m=l m=l 

where we have re-defined for brevity 

S2mi9) = '^ + S2m{g) . (3.11) 

4 Exact asymptotic expansion 

At this point, we start looking for the solution of equations (I3.10p for g -^ +00 in the 
form of an asymptotic series, i.e. 



°° (n) 00 (n) 

So™ , , . V ^ o-; 



S2r 



.(9)=E^. ■'— .w^E^' ("; 



n=0 ^ n=0 ^ 



Albeit we well know that the unique knowledge of the asymptotic expansion will never 
identify the function itself (contrary to the convergent weak-coupling series [9],[11]), we 
will see in the following many favourable issues, besides the contact with semi-classical 
string theory results. From the explicit solution of (I3.10p up to the order l/g"^, we have 
guessed the form at all orders: 
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with, implicitly, these expressions different from zero only if n < m. Moreover, the 
coefficients 6„ satisfy a recursive relation whose solution is more easily written in terms 

of their generating function 
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Consequently, we may also derive the explicit values 
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where E2k are Euler's numbers. To summarise equations (14.11 14.21 14.51) . we obtain the 
following strong coupling solution to the linear system (I3.10p : 



^^' ^ "" ^ ^ r(m - n + 1) ^2n-i Z^ (2A;)!(2n - 2A; - 1)! 
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As these series truncate, they are asymptotic in a very peculiar sense, namely they 
are correct up to non-expandable corrections which still go to zero (cf. below). From 
definition (13.71) and property (13. 4p we can read off 

s{0) = j^s,{g) = ^-Mg) , (4.7) 

which entails the exact asymptotic expression fi{g) = —1. 

It is very interesting to provide a possible summation for the expansions (14. 6 p 
(though finite!) by using an integral representation for the ratio of the two gamma 
functions and then by exchanging the sum over n with the integral, and finally sum- 
ming up: this procedure provides these new functions of g 
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Instead of the first generalised scaling function we obtain 

From this integral representation we see that the asymptotic behaviour becomes cor- 
rected by exponentially small terms, behaving as g^^'^exTp(—7ig/^/2). We will see in next 
section that this is not quite the right behaviour because of the wrong power g^^^. This 
is obviously connected to the wrong (convergent) weak coupling expansion of (14. 9p . 



5 Numerical solution: the non-asymptotic terms 

Remarkably, the linear system fl3.10p shows up the same matrix as in [15] for the BES 
equation, but a more involved forcing term. Explicitly, it may assume the form 

oo 

s,{g) = b,{g) - J^iK'^^ia) + 2Kl'i{g))sU9) , (5.10) 



m=l 



where 



Kl^^g) = 2iNZ),m, K^l{g) = A{PNZQNZ),^ (5.11^ 

h{g) = -{N + APNZQN)a^ - Na^ - A{NZP + APNZQNZP)^ , 



with 



A^ = diag(l,2,3,...), P = diag(0, 1,0,1...), Q = diag(l, 0, 1, . . . ), 

t = (l,0,l,0,...), a = (ai, 0,03,0,...), a = (0, 02, 0, 04, • • •)• (5-12) 

At least formally, the solution may be written as 

s{g) = (J+K('") + 2i^(^))-^6, (5.13) 



where I is the identity matrix. This similarity with ^15] is crucial for what concerns 
the numerical treatment, in that everything works well (as in [l5]), because (the matrix 
elements are the same and) the forcing elements have a similar behaviour for large g. 

In the following we will be interested in a numerical analysis of f l5.13p . emphasising 
the deviation of the strong coupling behaviour of the first few Sm{g) from the exact 
expansion of the previous Section. In this respect, a sufficient accuracy is achieved 
upon truncating at rra = 30: this restriction still yields reliable results up to about 
^~20. 

5.1 The strong coupling corrections and the string side 

Since the asymptotic expansion for any Sp{g) truncates significantly at the (inverse) 
power 1/g^, the additional (non-perturbative) exponential corrections are simply given 
by subtracting this to the numerical solution of fl5.13p [j. At any rate, our main concern 
here is Si{g) by virtue of its relation (14. 7p to the generalised scaling function fi{g), whose 
asymptotic expansion f{g) = —1 is in perfect agreement with the string energy density 
A — s = 7 + L (cf. (II. 2p ). as perturbatively expanded in [6l[H], because at present it does 
not seem to have the term L. Moreover, beyond this asymptotics there are additional, 
non-perturbative, exponential corrections, whose leading order (see below (15.140 ) is 
reported in figure [H where the large g, asymptotic behaviour is already reached for 
g ~ 3.5. 



^At the moment, we observe the same exponential behaviour, with different constant factors (cf. 
below), for the first few Sn{g), n = 1, . . . , 5. 

9 




Figure 1: Comparison between numerical solution (dots), fit as in formula (15.141) (dark 
grey line) and integral representation (14.91) (black line). Detail in the small window: 
comparison between numerical solution (dots), weak coupling expansion of [9l[ll] (light 
grey line) and integral representation (14. 9 p (black line). 



We can then proceed to a more quantitative analysis. At first, we perform a fit of 
the numerical evaluation of (15.131) . in order to check whether it is compatible with the 
proposed exact asymptotic expansion (14.61) . This analysis - performed at the moment 
on the first few Sn{g), n = 1, . . . ,5 - rules out the possibility that other terms with the 
form l/g'^ might appear. Then, we find reasonable to conjecture a non-perturbative, 
exponential term coming into play. In fact, a simple dimensional consideration has led 
the authors of [S] to claim that, regarding fi{g), the additional string energy density 
□ is simply proportional to (j and) the 0(6) sigma model mass gap, m, or in formulae 
{fi{g) + l)j = cm{g)j as long as j ^ m (still g -^ +cxd). |j As a consequence, we use 
the functional form implied by the mass gap formula ([8] and reference [13] therein). 



I.e. 



fiig) 



'l + ng^/'^e ^ 



g 



CXD 



(5.14) 



with only k as a free parameter to fit, and we find a perfect agreement with the data 



®This has to be summed to f{g) Ins and is due to the 5*0(6) charge j. 

^We may expect that m{g) in this formula differs from the mass gap as soon as j ^ g is no longer 
valid. 
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in the region g G [3,20], resulting in the best fit estimate k = 2.257 ± 0.009. On 
the other hand, the 0(6) small coupling calculations of [8] entail the exact value k = 
2^/^7r^/^c/r(5/4) = (2.265218666...) c, which compels us to the natural prediction that 

We can now comment a little about the integral summation (14.91) to fi{g)- Despite 
the fact that it is able to qualitatively capture the main features of the numerical solu- 
tion for a wide range of g, it still fails to reproduce both the weak-coupling expansion, 
and the power factor g* of the non-perturbative term. We have a clear picture of this 
situation in fig. [TJ Finally, we shall remark the perfect agreement between the numerical 
solution and the first terms of the weak-coupling expansion computed in [9l [TT] . 

6 Some final comments 

Unlike all the genuine (i.e. infinite) asymptotic series involved in the derivation of the 
leading term f{g) (see, for instance, [16]), which amounts to twice the cusp anomalous 
dimension [l7j, all the analogous series for fi{g) truncate, thus giving rise to polynomi- 
als in 1/g. Moreover, the latter become non-trivially corrected by exponentially small 
contributions. At this stage, a deeper comparison with string results ([B] and references 
therein) is really compelling, but still demanding because of the summation of all the 
logarithmic contributions in the string loop expansion [6]. Yet, this is the direct con- 
sequence of an ordering of different limits and the present results shed some light on 
what happens to the higher /„, n = 2, 3, . . . . 
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